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ABSTRACT 



This thesis contains a handout covering the Fundamental 
Physical Laws (Continuity, Momentum, and Energy) used in 
Aeronautical Engineering which are transformed from the con- 
trol mass or system form into the control volume form. It is 
intended that this handout serve as a self-studying guide for 
students in the core of the Aeronautical Engineering Program 
at the Naval Postgraduate School and as a reference during 
the graduate level courses. 
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I. INTRODUCTION 



A. THE PROBLEM UNDER CONSIDERATION 

In the Aeronautical Engineering Curricula at the Naval 
Postgraduate School the development of expressions for the 
Fundamental Physical Laws (Continuity, Momentum, and Energy) 
is covered at different depths in courses such as Basic Fluid 
Dynamics, Basic Thermodynamics and Gas Dynamics during the 
core . 

Due to the lack of time, these topics are usually treated 
in an introductory fashion and the notes taken by the students 
are necessarily incomplete with some repetition between the 
various courses. 

Also, in some cases the development of these equations is 
done without a systematical recognition that contro 1 -vo 1 ume 
approach is used in opposition to the system or fixed-mass 
approach (because of the inconvenience of the latter for prob- 
lems related to fluids). In addition to the problem previously 
described, it happens sometimes that the student is not aware 
that the equation the professor is using is a "disguised" form 
of one of the Fundamental Physical Laws for which he may al- 
ready have some intuition. 

The situation may be depicted in the following chart: 
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TOPIC: 


FUNDAMENTAL 
PHYSICAL LAWS 
INTEGRAL FORMS 


FUNDAMENTAL 
PHYSICAL LAWS 
DIFFERENTIAL 
FORMS 


PARTICULAR ,MS 

NAVIER-STOKES, 
ETC. 


DEVELOPED 

AND 

USED IN 
COURSES; 


2041, 2042 


3043 


3043, 4431 , 4432 
4501 , 4502 , 4511 
4512 


STAGE IN 
THE 

PROGRAM: 


CORE 


GRADUATE 


OBJECTIVES: 


CREDIBILITY 

^ AND ^ ^CONTINUITY OF SUBJECT MATTER-^ 

RIGOR 


< nMlMUil^L^j ’ 



B. TYPE OF STUDENTS ATTENDING THE CORE 

A considerable number of students taking the sequence of 
core courses are coming from the fleet after a number of years 
(between four and eight), since the completion of under- 
graduate education in Science or Engineering. 

It can be said that sufficient background exists in most 
of the cases but it is buried under several years of perfor- 
mance of nonacademic duties. 

The most noticeable effect of these inactive years is the 
difficulty in dealing with abstract concepts such as mathe- 
matical symbols and models intended to depict a situation of 
engineering interest (the control volume, for example). 

C. A PROPOSED SOLUTION FROM A STUDENT STANDPOINT 

The author believes that a handout received by the students 
before or during the core will provide them with a set of notes 
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for self-study during the core and for future re^^ererce in 
the graduate sequence of courses. 

The proposed handout on the Fundamental Physical Laws is 
a version of material covered in many books, some of them 
being more rigorous and elegant but often not appealing to 
the students. 

What is it that bothers a student in this situation when 
the topic is treated in class or in his textbooks? 

First of all, the terminology. The reader must be warned 
that the author is perhaps strongly biased in this respect 
because English is not his native language. But even students 
without this problem sometimes are confused by words or do not 
know the exact meaning of some terms in Thermodynamics or 
Fluid Dynamics. For example, look up the word "steady" in a 
dictionary. Reference 1 gives the following meanings for it: 
firm, stable, not shaky, regular, uniform, continuous, constant 
in behavior, etc. Assuming that the student knows all these 
meanings, if the word "steady" is used without rigorously 
defining it in a scientific context (i.e., partial derivatives 
with respect to time are zero), a sensation of ambiguity will 
remain and no convenient abstraction can be expected to last 
in the student's mind. 

Important terms must be carefully defined from the begin- 
ning and the other terms introduced as needed during the 
development. The use of these terms must be shown by means 
of examples and the degree of retention by the student must 
be frequently tested. (Notice how the term system is treated 
in the text . ) 
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A second source of trouble is the prolific contribution 
to Science by men such as Newton, Euler, Bernoulli (three of 
them! Jacob (1664-1705); Johan (1667-1748); and Daniel 
(1700-1782)); Laplace and many others. When someone quotes 
"Newton's Law," to which one is he referring — the First, 

Second, or Third? Maybe it is Newton's Law of Friction. 

Sometimes the student is taking, simultaneously, a course in 
Mathematics where the same names are also quoted. 

The fundamental points for a student in the core in this 
respect are: 

1. Not to get confused by the terminology and partic- 
ular names . 

2. To recognize where an equation originates; what is 
the significance of the various terms; and what are the con- 
straints and assumptions for its derivation. 

3. To make actual use of the equations. 

Another problem arises from Algebra. It is not that this 
branch of Mathematics is obscure or difficult for a student at 
this stage of his education, but that it is simply uninteresting. 

Oftentimes during a long and tedious derivation of an 
equation the student's mind wanders and the sensation of 
"magic" results when a final and neat expression appears. The 
following corrective actions may be proposed: 

1. State from the beginning the starting and ending 
points of the derivation. 

2. Underline the milestones in the derivation. 

3. Ask questions of the student in the middle of the 
derivation to keep him interested or simply awake. 
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4. Skip on purpose some steps but ask the students 
to go through them for completeness. 

5. Use the most compact notation available but only 
after being sure that the notation is properly introduced and 
that this greater degree of abstraction is consistent with the 
level of the student. 

Perhaps many of these problems are avoided by using a 
textbook writing technique called "Programmed Instruction." 

The opinion of the author is that writing a handout fully 
using this technique is out of the scope of a Master's Degree 
Thesis, but some of the features of this technique can be 
incorporated and the possibilities of improvement in this 
respect are unlimited. 

D. TUE VARIABLE OBJECTIVES OF THE HANDOUT 

The objective of the proposed handout is not necessarily 
to teach a technique or a skill (such as the use of the Laplace 
Transform to solve certain differential equations). The objec- 
tives of this handout are variable from the student standpoint, 
accordingly to what stage of his education he is at. 

At the beginning, during the core, the main text of the 
handout is intended to satisfy the objective of CREDIBILITY 
and RIGOR. The student should be familiar (and if he has for- 
gotten, the very first part of the handout will serve as a 
refresher) with the system or fixed-mass approach. The hand- 
out will lead him neatly and rigorously to the control volume 
approach through a transformation or extension of the system 
approach . 
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This needs to be done only once and it can be done using 
a generalized notation. Also it is believed that a self- 
study text has the poten ti al of removing the usual class- 
room problems. 

In the next steps the student is led through specializa- 
tions of the control volume formulation. The students who 
are in the core are not expected to carry out these further 
specializations, but he is expected to do so in his graduate 
courses. Here, the objective of the handout is to provide 
CONTINUITY OF SUBJECT MATTER through the whole sequence of 
courses . 

The same objective is pursued in the extensions leading 
to the differential form of the Fundamental Physical Laws and 
to other particular forms. He is not expected to carry out ^ 
these derivations, but he is expected to understand where the 
equations come from and (eventually) the use of them in the 
context of the graduate level courses. Here the "mystery" of 
the equation written down without proof might be removed and 
the student uneasiness with an advanced subject alleviated. 

The entire handout should serve properly as a review and 
reference during all the stages of his career at the Naval 
Postgraduate School. The objective in this respect is HANDI- 
NESS. This is perhaps the most ambitious aim of this work. 

E. THE LEARNING PROCESS 

This thesis has been written according to the following 
basic ideas with respect to the learning process; 

1. The assumption that the student is motivated; he has 
the desire to learn. 
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2. No mathematical background is assumed beyond undergraduate 
level . 

3. Certain basic vocabulary is assumed but the important 
concepts are defined or explained carefully and their 
understanding emphasized and tested in several parts of 
the text. 

4. There would be systematic reading of the context during 
the self-study use of these notes. The student will not 
proceed to another section without a fair understand ing of 
the previous. (Approximate relative times for reading 
are stated . ) 

Along the text, several quizzes, programmed questions and 
exercises are inserted to help the student recognize his own 
degree of retention of what he was supposed to learn. 

When a nevj quantity is introduced and assigned a symbol, 
the correspond! ng units are shown in the English system only 
as an example to help the reader recognize and associate the 
nature of the quantity. 

Also, to some restricted extent, some provisions are made 
in order to allow the student to extend his recently obtained 
knowledge to new situations. 

F. SUMMARY OF THE PRESENTATION OF THE SUBJECT IN THE HANDOUT 
The general idea followed during the development of the 
thesis was to treat the topic at the highest degree of 

generality and compactness consistent with the level of the 

< 

student and then to particularize to arrive at useful forms 
of the equati ons . 
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It can be mentioned in this respect that the reduction 
of the expressions of the three Fundamental Physical Laws to 
one using a generalized notation, was based on Ref. 2 after 
consulting several undergraduate level books and notes 
[Refs. 3 through 11]. 

The concepts involved in the control volume approach were 
applied to this generalized expression of the Fundamental Phys- 
ical Laws for a system. Later a general expression of the 
Fundamental Physical Laws for a control volume was developed. 
Then this expression was specialized to the particular laws in 
integral form and later in differential form. 

The advantage of this method is that the mathematical 
transformation needs only to be done once, enabling the student 

to concentrate on the actual mathematics without having to 

0 

identify the "proof" with a single Fundamental Physical Law. 

The steps in particularization of the equations are con- 
tained in the supplements where the Navier-Stokes, Bernoulli 
and Euler equations are developed having in mind that very 
often a fair knowledge of these equations is required in 
courses such as Boundary Layer Theory, Convective Heat and 
Mass Transfer, Magnetohydrodynamics, Advanced Gas Dynamics, 
Turbomachinery and Aerodynamics of Wings and Bodies. 

G. THE SCOPE OF THE APPROACH CHOSEN 

This thesis is a student's point of view on the subject, 
taking into consideration the difficulties found during his 
education . 



11 



This type of approach may be helpful to fellow studer^s 
in the same v/ay as it has been for the author. 

Under the title of this thesis, a broader range of sub- 
jects can be covered. However, some amount of limitation was 
found to be necessary in order to fulfill most of the purposes 
enunciated earlier in this Introduction. 

It is necessary to mention that the treatment of the 
Equation of State which is complementary to the Fundamental 
Physical Laws has been omitted. The reason for this is that 
this topic is covered in a thesis written at the same time 
by another student [Ref. 12]. 

II . THE HANDOUT 

In order to have the proposed handout in a style and 
format more useful for its purpose, it is inserted in the 
Appendix. 

Ill . A SURVEY OF THE USE OF THE HANDOUT BY STUDENTS 

An earlier version of the handout was distributed to test 
its acceptance by the students. 

This was done in Course AE 2041, Basic Thermodynamics, 
where Ref. 6 was used as a textbook. Upon the completion of 
the Course the instructor made a survey of the opinion of the 
students about the handout. The following are some results 
of such a poll: 
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95 % read the handout in its entirety. 

88 % had sufficient background to understand the 
material . 

95 % felt that the material was of sufficient 
difficulty to merit a handout, 

95 % felt that the material was presented in a 
logical order. 

When asked whether they would prefer their textbook, the 
handout, neither, or both, the response was: 

40 % present text 
35 % handout 
15 % both 
10 % neither 

One important conclusion based on the survey is that at 
least the handout received was a desired instructional aid. 

With the improvements made since then and other improve- 
ments that can be made after testing it again, the handout may 
be transformed in a very useful instructional aid to be re- 
ceived by the students at the beginning of their career at 
the school. 
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IV. RECOMMENDATIONS FOR FUTURE WORK 



The handout may be improved mainly in the supplements. 

In Supplement C, the differential form of the Energy Equa- 
tion is only stated, suggesting that its development is simi- 
lar to the development of the Navier-Stokes Equations which 
in turn are not derived in detail for the sake of conciseness. 

Perhaps a few more steps may be added in the development 
of the Navier-Stokes Equations. The complete derivation of 
the Energy Equation in differential form is needed to satisfy 
the requirement of completeness for reference use. 

There are also opportunities for improvement in the main 
body in the style of the presentation. After being used by 

o 

the students, the programmed characteristics may be increased 
to overcome difficult steps. 



14 




# 



APPENDIX 



Naval Postgraduate School 
Department of Aeronautics 
Monterey, California 



THE CONTROL VOLUME CONCEPT 
IN AERONAUTICAL ENGINEERING 



From an M.S. Ae . E. Thesis 

By: LT Carlos Tromben 

Thes i s *Advi sor : 

Asst Professor Oscar Biblarz 
1971 



15 



INDEX 



I. PRELIMINARIES - 20 

A. A SUGGESTION ON HOW TO USE THIS HANDOUT 20 

1. Self-Check Quizzes 20 

2. Programmed Questions 20 

3. Exercises 20 

B. FUNDAMENTAL PHYSICAL LAWS 21 

C. BASIC DEFINITIONS 21 

1. A word about the Notation 21 

2. A word about FUPLAs 23 

D. SELF-CHECK QUIZ NO. 1 26 

II. DEVELOPMENT OF EXPRESSIONS FOR THE FUNDAMENTAL 

PHYSICAL LAWS FOR A SYSTEM 27 

0 

A. CONSERVATION OF MASS --- 27 

B. MOMENTUM EQUATION 28 

C. ENERGY EQUATION 30 

D. SUMMARY - 32 

E. SELF-CHECK QUIZ NO. 2 - 35 

III. DEVELOPMENT OF EXPRESSIONS FOR THE FUNDAMENTAL 

PHYSICAL LAWS USING THE CONTROL VOLUME APPROACH - 37 

A. THE DEVELOPMENT IN GENERALIZED NOTATION 37 

B. PARTICULARIZATION OF THE GENERAL EXPRESSION --47 

1. Continuity Equation 47 

2. Momentum Equation 47 

3. Energy Equation 47 

C. SELF-CHECK QUIZ NO. 3 49 



16 






I 



« 



IV. 



EXTENSION OF THE CONTROL VOLUME APPROACH TO 
OTHER FUNDAMENTAL PHYSICAL LAWS 51 

A. ANGULAR MOMENTUM 51 

B. SECOND LAW OF THERMODYNAMICS 52 

C. SELF-CHECK QUIZ NO. 4 54 

V. APPLICATION TO PARTICULAR PROBLEMS 56 

A. CONTINUITY EQUATION 56 

B. EXAMPLE OF APPLICATION OF THE SYSTEM AND 
CONTROL VOLUME APPROACH TO AN UNSTEADY 

PROBLEM 58 

1. System Solution 

2. Control Volume Solution 60 

VI. SUPPLEMENT A: BODY AND VOLUME FORCES 64 

VII. SUPPLEMENT B: INTERMEDIATE STEPS IN THE 

DERIVATION OF THE ENERGY EQUATION 66 

VIII. SUPPLEMENT C: EXPRESSIONS FOR THE FUNDAMENTAL 

PHYSICAL LAWS (FUPLAs) IN A DIFFERENTIAL FORM -- '' 

A. CONTINUITY - 72 

B. MOMENTUM 73 

C. ENERGY - 78 

D. SELF-CHECK QUIZ NO. 5 - 80 

BIBLIOGRAPHY - - 82 



17 



I 






TABLE OF SYMBOLS AND ABBREVIATIONS 



A area on the control surface 

CV control volume 

/N /N 

CS control surface 

E energy 

e energy per unit mass 

g acceleration of gravity 

g dimensional constant in Newton's 2nd Law 

^c 

h enthalpy per unit mass 

h^ total enthalpy per unit mass 

k thermal conductivity 

m mass 

moment 

n unit normal vector 

^ momentum 

P rate of work 

p pressure 

Q heat 

Q heat flow. Heat per unit time. 

q heat flow per unit area. Heat per unit time per unit 

area . 

r position vector of the center of mass of a system of 
particles or position vector of a particle. 

S entropy 

s entropy per unit mass 

syst system 

T temperature 

t time 
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U internal energy 

V volume 

V velocity. (v = |v| speed) 

W work 

X an extensive property that is "conserved" within the 

system 

V a quantity that represents how the system interacts 
with the surroundings 

z height above a datum point 

u viscosity 

C arbitrary vector 

p density, mass per unit volume 

$ dissipation function 
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I. PRELIMINARIES 



A. A SUGGESTION ON HOW TO USE THIS HANDOUT 

This handout is supposed to promote a better understanding 
of the Fundamental Physical Laws (Continuity, Momentum, and En 

ergy) and is to be read out of the classroom (i.e., self-study) 

To obtain better results, (usually) at the beginning of 

each chapter or section there is a note suggesting how far 
you, the student, should go on reading continuously. We do 
not forecast how long it will take to read and understand each 
unit because of individual differences, but our guess is that 
the first two units will take you about one hour. The esti- 
mates on the following sections are indicated relative to 
the first. 

If you do not have the time now to read continuously as 
suggested, try it some other time, because partial reading 
will not benefit you in any way. 

Also, there are th ree types of testing material inserted 
along the text and you should understand their purposes: 

1 . Self-Check Quizzes 

The purpose of these quizzes is to help you recognize 
your own degree of retention before proceeding to the next 
page of the material. Be sure to correct the wrong answers, 
and do not continue unless you feel you understand the dif- 
ficulties you had. The answer is usually at the bottom of 
the page, written upside down. 

2 . Programmed Questions 

The primary purpose of these questions is to remind 
you of important points presented in previous sections. You 

must fill the blank space left on purpose and the answer is 
usually at the bottom of the page, written upside down. 
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3 . Exercises 



These are not truly testing material. Their r-urpose 
is to make you think about further implications of the subject 
in an open-ended fashion. No ansv/er is given for these. 

B. FUNDAMENTAL PHYSICAL LAWS 

In the solution of problems in Aeronautical Engineering, 
concepts from the Physical Sciences called "Fundamental Phys- 
ical Laws," must be applied. These laws relate to mass, 
momentum, and energy and are commonly referred to as "conser- 
vation laws" because under many circumstances they are truly 
conserved; however, the word conservation in its common usage 
takes a broader meaning which may be stated as "accounted for" 
or "not destroyed." We will use the term "Fundamental Physical 
Laws (FUPLAs)" because this will permit us to include the second 
law of thermodynamics relating to the entropy which is not con- 
served in a natural process. In addition to the FUPLAs, in- 
formation is needed about the fluid itself usually contained in 
the Equation of State. This equation describes a set of unique 
relationships between state properties and may be found in 
algebraic form (the ideal gas law), tabular form (the steam 
tables), or graphical form (p-v-T surfaces, etc.). 

C. BASIC DEFINITIONS 

Read continuously from here to the end of Section D. 

1 . A Word about the Notation 

Through all this work a (caret) will be used over 
the first letter of a word which has a specific meaning in 
this text. For example, whenever the word cycle appears with 
a '' over the letter "c" (i.e., cycle), means a precise concept 
and that there is no other meaning for the same word. 

All words beginning with a letter accented with ^ are 
carefully defined in the following list. 
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Fluid; a substance that deforms continuously under the action 

of a shear stress. 

/\ 

System; a collection of matter of fixed identify (Control 
Mass ) . 

Control Volume (CV): a region in space of fixed size through 

which fluid flows. It may be stationary, moving uniformly 
(inertial), or accelerating ( non- i nerti a 1 ) . In this work we 
will consider the CV always fixed with respect to an inertial 
coordinate axis system,^ unless otherwise specified. This 
volume is bounded by a control surface. 

We want to insist on the idea that, although the con- 
trol volume has to be fixed with respect to an inertial sys- 
tem of coordinates, it does not necessarily need to be fixed 
in space. The following illustration shows this. The vehicle 
may be moving through space but the control volume is fixed 
with respect to the vehicle. 




Figure 1. Fixed Control Volume with respect to inertial 
frame of reference. 



1 

Our development is constrained to this case for simplic- 
ity and to fit into the scope of the handout like this. Impor- 
tant Aeronautical Engineering applications require non-inertial 
frames of reference, such as in turbomachinery, for example, 
[Ref, 13]. 
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/N /s 

Surroundings: everything outside the system 

^ /S 

Extensive Property: a property that is mass dependent. It 

can be a scalar, as energy or a vector, as momentum 

Fundamental Physical Laws (FUPLAs); there are three funda- 
mental physical laws which, with the exception of rela- 
tivistic and nuclear phenomena, apply to each and every flow, 

independently of the nature of the fluid under consideration 

^ 

State: the condition of the system characterized by the 

values of its properties 

Process: the path of the succession of states through which 

the system passes 

Cycle; a system that undergoes a series of processes and 
always returns to its initial state is said to have gone 
through a cycle 

Properties: characteristics of a system that define its state. 
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TABLE I 



FUPLAs (From Observation) 


MATHEMATICAL FORMULATION 


Conservation of Mass 
Newton's Second Law of Motion 
First Law of Thermodynamics 


Continuity Equation 
Momentum Equation 
Energy Equation 



Note that Table I is incomplete. Other FUPLAs will be seen 
later. 

2 . A Word about FUPLAs 

Perhaps the reader has become uneasy about the previous 
wording. All those terms, with the caret (which denotes a 

/N /N 

specific meaning) and the solemn title "Fundamental Physical 
Laws," may lead one to think that those laws have something 
mysterious or difficult to understand. 

This should not be so. 

These laws were deduced from observations of nature. 
Nature behaves this way. It has always behaved so and nobody 
has reliably observed a contradiction to this behavior. 

This is the only "proof" of these laws. 

/S /N /\ 

Later these Fundamental Physical Laws are put into a 
mathematical form in order to manipulate and work with these 
concepts in the precise and synthetic language of the mathe- 
mati cal symbol s . 

This almost trivial example shows that the FUPLAs 
belong to everyday experience: If you fad a dozen apples 

which you were to squeeze to make apple sauce with, you would 
not be conserving the number of apples. What is to be 
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accounted for is the matter (mass) when the seeds and other 
discards are taken into account. 

FUPLAs are the generalization of many reliable observa- 
tions (i.e., matter is not destroyed, it is transformed into 
other forms of matter or into energy). 
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D. SELF-CHECK QUIZ NO. 1 



Select the proper answer 

1. The .Control Volume is: 

a. a collection of mass of fixed identity 

b. always the volume occupied by the system 

c. a region in space of fixed size through which 
fluid flows 

d. always stationary 

/N /X /\ 

2. The Fundamental Physical Laws: 

a. are deduced from observation of nature 

b. can the proven using highly complicated mathemat- 
ical procedures and as a result of this, most of 
the textbooks do not have them 

c. apply to certain types of fluids, flowing under 
speci a 1 condi ti ons 

d. have no mathematical expressions 



(e )2 pue (D)i :sj3Msuv 



EXERCISES: 

What is conserved when an object of mass m falls from a 
height h in a gravity field? 

A typical reference state for air is standard conditions. 
Do you know what this means? 
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1 1 . DEVELOPMENT OF EXPRESSIONS FOR THE FUN DA'iEriT^-L PHYSICAL 

LAWS FOR A SYSTEM 



Read continuously up to the end of Section E. 

A. CONSERVATION OF MASS 

A system is 

Let us consider the following system 




Figure 2, Gas in a closed cylinder 

Experience tells us that matter will not flow in or out 
through the solid boundaries of this system (or any other), 
as time elapses or as the identified system wanders through 
space. Also, matter will not be created spontaneously within 
the system. From' this last statement we may write 



msyst ” constant 

m.„_. = / .dm 

syst syst 

paxig go sseiu go /Cgiguenb b 



( 1 ) 

( 2 )’ 

:U3MSU\/ 
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If p = ^ , then - ^"syst“‘^'' 

Differentiating v/i th respect to time in Eq. 3 and Keec'ng 
in mind that 'the left-hand side is a constant (from Eq. 1) 
we get ^ 



dt Jsyst 



pdV 



(4) 



Note that in the above, as well as throughout the rest of 
this Section II, integrals over the mass (as in Eq. 2) are 
being used. They are changed into volume integrals provided 
we can use the mass density p. In other words, for the system 
analysis, integrals over a fixed quantity of mass are being 
dealt with, and eventually these basic laws in terms of control 
volume will be expressed. In order to avoid confusion as to 
which one the volume integral belongs, the term "syst" will 
appear by the integral sign when a fixed mass is being analyzed 
and the term "CV" when a fixed or control volume is being 
analyzed. 



B. MOMENTUM EQUATION 

Consider a differential volume element of our system which 
is traveling at velocity v. This volume element will have a 
mass dm and again 



dm 


= 


pdV 






(5) 


and a momentum 












d^ 




vdm 






(6) 


dP 


= 


vpdV 






(7) 


So that the 


total momentum 


for the system will be 






r 






r 






syst 


d^ = 


syst 


vpdV 


(8) 






Differentiating Eq. 8 with respect to time 
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r 



dt 



syst 



dP 



dt 



J 



syst 



vpdV 



(9) 



The Momentum equation is the expression of Nev^/ton's Second 
Law of Motion which states: 

/s 

"The time rate of change of momentum of a system 
is proportional to the net external force acting on 
system and takes place in the direction of the net 
force . " 



(Is there a "proof" for this law?) 
Therefore, ^ 

F 



net 



dt _Jsyst 



dP 



(10) 



The net force is made up of contributions of two kinds of force, 
namely. 






F 



net 




vol ume 




surface 



( 11 ) 



The above distinction regarding the net force is needed 
because of the different nature between, say, the force of 
gravity and a pressure force (Supplement A expands on the nature 
of these forces ) . 



Inserting Eq, 9 into Eq, 10: 



r 



net 



dt 



J 



syst 



pv dV 



If there is no net external force, Eq. 8 becomes 

r 

d 



dt 



sys t 



dP 



( 12 ) 



(13) 



or 



P = constant (i.e., momentum is conserved) 

(14) 



The student should eventually learn to recognize the dif- 
ference between the collision of two bill.iard balls and of two 
balls of putty in terms of the momentum equation. 
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Again, all that has been considered so far, with respect 
to the Conservation of Momentum, is restricted to the syste. 
approach. At this point the above will not be very useful 
if the system is difficult to identify. 



C. ENERGY EQUATION 

The Energy Equation is the mathematical expression of the 
First Law of Thermodynamics, which can be stated as follows; 

"If a system is carried through a cycle, the total heat added 
to the system from its surroundings is proportional to the 

^ /N 

work done by the system on its s urroundings . " 

Let us have a system, that is, a quantity of matter of fixed 
identity, traveling through space (changing its position with 
respect to an inertial frame of reference as time elapses). 

/N 

Let the path of the system be in such a way that it will re- 

turn to its initial state. In doing so, it has performed a 
/\ 

cy cl e . 

/S /N 

During the cycle, heat may be added to the system by its 
surroundings and in turn, work may be done by the system on the 
surroundings . A useful cycle is depicted in Figure 3. 



Notice that the previous discussion and Figure 3 refers to 

a system undergoing a cycle, that is, 

. A similar discussion may be made for a 

system undergoing a single process. 



Since energy is related to heat and work, for a process — 
an expression for the First Law of Thermodynamics in the 
rate form can be written; 



Rate of heat 
flow into the 
system 




rate of work 
done by system 
to the 

s urroundi n gs 




rate of change of 
energy within the 
, system 



(15) 



Here, the following sign convention is used: Heat added to 
the system by the s urroundi n gs , positive (+). Work done on 
the surroundings by the system, positive (+) . 



•eqeqs [euLbLuo sqi oq luaqs^s eqq 
J .0 aqeqs eqq sujnqou s/CeM[? qeqq sassaoojd go saiuas e ;jaMSuv 
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Figure 3. Gas cycle 




WORK 

W ( + ) 



Figure 4. Sign convention for the First Law 
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Some authors use 



This sign convention is not universal, 
other conventi ons . 

Let us define the following symbols: 

Q rate of heat flow (e.g., in 

IP rate of work done (e.g., power i.e., in or watts) 

e total energy per unit mass (e.g., 

Then, one may write the energy equation for a system from Eq. 
15; r 



• 

Q 


- JP = 


d 

dt_^ 


syst " 


(16) 






r 




• 

Q 


- IP = 


dt 


syst ^ 0 '1'' 


(17) 






The idea that energy cannot be created or destroyed within 
a system, but that it can be exchanged is an old one and is 
easily accepted due to the fact that it has never been ob- 
served to be violated in nature. 

The problem again is to extend the Energy Equation to cases 
where it is difficult to follow the control mass. 

D. SUMMARY 

The expressions for the FUPLAs for a system is summarized 
in Tables II and III. 

Table II will be extended to other Laws which has not yet 
been discussed. 

In Table III there are some symbols not defined so far. 
Refer to "Extension of the Control Volume Approach to other 

/N /N /N 

Fundamental Physical Laws" for explanation (Chapter IV). 
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TABLE II 



Fundamental . 
Phys i cal 
Law 


Name Usual ly Gi ven 
to the Equati on 


Expression for 

/N 

a System 


Conservat i on 
of Mass 


Continuity Equation 


d 

dtj 


syst P'*'' = » 


Newton ' s 
Second Law 
of Motion 


Momentum Equation 
or 

Momentum Theorem 


dt 


r 

. vpdV = F . 
syst ^ net 


First Law 
of 

The rmodynami cs 


Energy Equation 


d 


r 

syst 'P'*'' ' ^ - P 



TABLE III 



/s 

Fundamental 

0 

Phys i cal 

/N 

Law 


Name Usually Given 
to the Equation 


Expression for 
a System 


An gul ar 
Momentum 


Moment of Momentum 


dt 


r 

5ySt(P«)P<iV , 


Second Law 
of 

Thermodynami cs 


Entropy 


dt 

J 


r r 

• 

syst^ syst? " 

J 



Looking at the two previous tables, a clear pattern is seen 
in the form of the expressions for the FUPLAs. We may gener- 
alize this form as r 



dt 



syst 



X dV 



Y 



X represents an extensive property that is "conserved" 
the system provided that the right-hand side of Eq. 18 
zero (Y=0 ) . 



(18) 

within 
i s 



Y is a quantity that represents h o^ the system interacts 
with the surrounding. Y's act on the system to change the X's. 

Note that there has been no difference in our description 
for a system composed of a differential quantity of 
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fluid or, say, a solid such as a cannon ball. We find t e 
system description, however, inconvenient in a flow situation 
because it is difficult to visualize it. 

Suppose you can identify a fixed quantity of air flowing 
into a fan jet aircraft engine. It would be very difficult to 
follow that identified mass of air as it travels through the 
engine and experiences changes in properties. See Figure 5 
below. 




Figure 5. A jet engine 



It would be much easier to choose a convenient control 

/N 

volume, fixed with respect to the engine for example and then 
investigate the changes in the air properties as it flows in 
and oiit of the control volume. 

Exercise: 

What would be a convenient control volume to describe what, 
happens to air as it passes through a smoking pipe? 
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E. SELF-CHECK QUIZ NO. 2 



Select the proper answer 



1 . A system is: 

a. the quantity of matter inside the control volume. 

b. a collection of matter of fixed identity. 

c. a collection of mass of fixed position. 

d. the set of substances that makes up a mixture. 

2. The Momentum Equation is: 

a. the expression of Newton's Third Law of Motion. 

b. the mathematical expression of a Law based on 
observation of nature. 

c. other name given to the Moment of Momentum Theorem. 

d. the expression of the Law of Conservation of 
Angular Momentum. 

3. The First Law of Thermodynamics; 



a. cannot be applied to a flow situation. 

b. has a mathematical expression usually called the 
Entropy Equation. 

c. is applicable to the interchange of heat and work 
between the system and the control volume. 

d. has a mathematical expression called the Energy 
Equation. 

/\. /N /N 

4. The expressions for the Fundamental Physical Laws for 
a system: 



a. show a clear pattern that makes possible to 
generalize them using a simple notation. 

b. are difficult to apply in certain flow situations. 

c. involves vector and/or scalar quantities. 

d. all the above are correct. 



•(p)fr ‘(p)e ‘(q)2 ‘(q)i -.suaMsuv 
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In the generalized expression for the Fundamental 
Physical Laws for the system the element of volutie 
dV is rel ated to : 

a. the Control Volume. 

b. a volume inside which properties change from 
point to point. 

c. the space occupied by the system which is 
stationary. 

d. none of the above is correct. 



(P)S 
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III. 



/N 



DEVELOPMENT OF EXPRESSIONS FOR THE FUN DA:-"'iT AL 
LAWS USING THE CONTROL VOLUME APPROACH 



F- ''Si CAL 



Read continuously up to the end of Section B, but only if 
you are acquainted with the material up to here. Otherwise, 
do not proceed until understanding the system approach. You 
can expect the reading of this unit to take about twice as 
long as each of the previous units. 



A. THE DEVELOPMENT IN GENERALIZED NOTATION 

Instead of developing separate expressions for every 

^ /\ /s 

Fundamental Physical Law, all of them can be done at once, 
using the generalized notation 



r 



dt 

J 



syst 



X dV = Y 



(18) 



Remember that in the above expression: 

1. Is valid for a only. 

2. XdV is an extensive property--a dependent quantity 

that is conserved the system when Y is zero 

(except for the entropy). 

3. Y is a quantity that affects the state of the system. 

4. The differential dV refers to the volume filled by the 

DO NOT CONFUSE WITH THE 

RIGOROUSLY AND COMPLETELY DIFFERENT CONTROL VOLUME. 

We will use a control volume (CV) fixed with respect to 

. Fluid can flow through this 

control volume. A typical CV is depicted in Fig. 6. 



iua:^sX’s aqeuipjooo [eiquaui 
iua:isX's uo sseiu louquoo 

N/ V ^ 

ULLi:}LM 

sseuj 
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Figure 6. A Nozzle as an Example of Control Volume 
I't is possible to proceed vnth our development using 

/V /V 

the above shown control volume, but a more general form or 
shape may be preferred, since we are trying to be general 
in the notation. In solving particular problems using the 

y\ /V 

control volume approach a wise choice of a convenient 
control volume, suitable to particular situations, has to 
be made. 
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Let 


u s 


con s i de r 


the f 0 1 1 ov/i n g 


general 


f 1 ovj : 


At 


t i me 


"t" the 


cont rol vol ume 


is e n t i 


irely occj 


sys 


tern 




Vcv) • 







d by 




Figure 7. Snapshot at time "t" 



39 



At time "t + ut" we can distinguish three different 
definite regions: 




»I>W 

r. s--' 

■ - 

v:c>. V 



Figure 8. Snapshot at time "t + At" 



Note that the system may have changed shape at time t + At. 
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Region I: common to the control volume and the system 

Region II; occupied by the portion of the system that left 

/V 

the control volume during the interval At 

/N /V 

Region III: portion of control volume unoccupied by matter 

belonging to the system because it changed 
position during the interval At 

/\ 

From calculus, the time rate of change of the extensive 

/N /V 

property X in the system is 



dt 






syst'^ = 1^'^t-O Ati 




syst 

@t+At 



X dV 




X dV 



(19) 



Remember that in this case the differential dV refers to 
the volume occupied by the system not the control volume 

Now at time t + At 

syst = region II + region I (20) 



0 r 



and at time t 



syst 



syst 



region II + CV - region III 



CV 



( 21 ) 

( 22 ) 



Using these last volumetric considerations; 



dt 






syst 



X dV=l i m 



At-^0 




r 



RI I 
t+A t 



XdV+ 






CV 

@ t+At 



XdV- 



RIII 
t+At 



XdV- 



CV 

@t 



XdV 



(23) 



Rearranging: 

r 



dt 



syst 



X dV= 1 i m 



At^O 




r 



X dV- 



0t+At 



r 






CV 

0t 



X dV + 



r 



RI I 
>>0t+At 



X dV- 



RIII 
0 t+At 



X dV 



(24) 
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Let us v;ork with the right-hand side of Eq, 24 taking two 
terms at a time. 



The first two terms represent changes with respect to time 
of the extensive property X inside the control volume. 



Thus using the definition of Ijmit: 

V 



1 i m 1_ 
At^O At 




r 



X dV - 



X dV = ^ 



3t 



X dV 



CV 



(25) 



/CV / 

Note that fhe right-hand side of Eq. 25 has a partial 
derivative, meaning that there is a change with respect to 
time of a quantity that depends on several independent vari 

/N. /N 

ables. In the control volume X = X(x,y,z,t). 



We will take care now of the last two terms of Eq. 24 



lim 1_ 
At-i'O At 



RII 

@ t + A t 



X dV 



r 






RII I 
IP t + A t 



X dV 



These terms, associated with regions II and III, represent 



the amount of the extensive property X inside these regions 
at time t+At . 



The amount of the extensive property X inside II at time 
t+At can be related to the amount of X that can enter during 
At to this region, crossing that part of the control surface 
represented by segment 3-4 in Fig. 9a. 

Similarly, the amount of property X inside III at time 
t+At can be related to the amount of X that can enter to this 

/N 

region during the interval At crossing that part of the control 
surface represented by segment 1-2 in Fig. 9a. 



Later we will find what this relation is, in order to 
transform the two last terms in Eq. 24. 



Referring now to Fig. 9b: Identify the system drawing a 

segmented line as in Fig. 9b. Also identify the segments 1-2 
and 3-4 of Fig. 9a. in the blank circles in Fig. 9b. (Answer 
in Fig. B2) 
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Control Vol’jme 
/\ 

System 



Figure 9a. One dimensional version of Figure 9b 







Figure 9b. Two dimensional projection 
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Let us find the mathematical expression of the 
property X flowing in and out of the control volu ;e ccruss the 
control surface. If dA is an increment in the control surface 
represented by segments 1-2 and 3-4 in Fig. 9, then 



V • n is the component of v perpendicular to dA, 

(v • n)dA is the volumetric flow rate and 
X(v • n)dA is the flow of property X 

Note that X(v • n)dAAt = amount of X that crosses dA in time At 

Since dV = dA(v • n)At 



Replacing this last expression for dV in the last term of 



the right-hand side of Eq . 24 



1 i m 1_ 
t-^0 At 



r 



RII 
@t+At 



X dV 



r 



1 i m 
t->0 



1 _ 

At 



X(v 



RII 
0t+At 



n ) dA 



At 



The integral in the right-hand side of the above equation 
represents the extensive property X that flowed into region II 
during the interval At and it was said before that this flow is 



across that part of the control surface represented by 3-4 in 
Fig. 9. Then defining the area corresponding to this 

segment, one ma}/ write 



lim 1_ 
t->0 At 






RII 
0 t+A t 



X(v • n) dAAt = 



A 



X(v • n) dA 
out 



(26) 



Now if a vector n' is defined as inward normal (inward and 
outward refer to the control surface) the quantity 



AinX(v • n')dAAt will represent the amount of property X 



that has flowed into region III during At, so 

r r 



1 im 1_ 
At-^0 At 



RII I 
0t+At 



X dV = 



Nn 



X(v • n' )dA 



and from Eqs. 26and 27 



lim 1 
At->0 At 




r 



RII 
0 t + At 



X dV - 



RIII 
M t+At 

r 



X dV 



(27) 



r 






A 



X ( v*'n ) dA 



out 



J 



^•n 



X(v • n')dA 



.( 28 ) 
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but n = -n' (by construction) 



lim ]_ 
At-^0 At 




dV - 



r 

RI II' 
@t + A t 



r 



-^^out 



X(v 

r 

^ ^in 



n ) dA 



X(v • n)dA 



(29) 



The combination of the two terms on the right-hand side 
of Eqs. 28 and 29 represents the flow across the boundaries 
of the control volume that occupies regions III and I at 
time t + At. We can rewrite Eq. 29 as 



lim 1_ 
At^O At 




dV - 



J 



Ri I r 
@t + At 




X(v • n)dA 



(30) 



Because the integration over Ain and Aout is the same as the 
integral over the whole control surface and nds e dA we can 
substitute Eq. 30 into Eq. 25 











r , 


r 




d 

dt 


^X dV 
syst 


_ 9 
9t 


cvX 






J 




J 


J 




which can 


also be 


wri tte n 


as 








r 


r 


r 



dt 




X dV 



J 



AW + 

CV 9t 




(31) 



(32) 



(Through the use of Leibnitz rule which is applied to the 

/N /N 

control volume in space.) 

The first term of the right-hand side may be considered 
as the rate of change of the quantity of X stored in the 
control volume while the other term as the net flux of X out 
of the control surface. 
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Equations 31 and 32 are very important because they re- 
late the time rate of change of property X within a systen 

/V 

or control mass with changes happening in the control volune. 



Let us examine the first term of the right-hand side of 
Eq. 31. It represents the time rate of change of the exten- 
sive property X within the control volume. It is a triple 
integral or a volume integral, so when we write 



r 



rn 



9t 






cv 



X dV 



we mean 



9t 



CV 



X dV 



This term is often called "local time rate of change" or storage 
te rm . 

The second term of the right-hand side of Eq. 31 represents 
the net efflux (outflow) of the extensive property X across the 
boundaries of the control volume, also called control surface. 

->• /s 

The vector dA is perpendicular to the control surface at every 
point of it and is positive when directed outward. The dot 
product V times dA gives the projection of v over the normal 
to the surface direction. Again: 

r rr 




This term is often called "the convective term." 

Equations 31 and 32 are transformation equations that en- 
able us to obtain an expression for the FUPLAs for a control 
volume from the corresponding equation for a system. 

From Eqs. 31 and 18 we obtain 




(33) 



This is the generalized equation for the FUPLAs in 

/N /N 

control volume form. 
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B. PARTICULARIZATION OF THE GENERAL EXPRESSION 



Now that Eq. 33 gives a general expression for the 
FUPLAs in terms of X and Y, it may be applied to specific 
cases: 



Continuity Equation 
X = p 

Y = 0 r 

3 



* SOS'" “ ■ 



(34) 



or 




2 . Momentum Equation 

Note that here both X and Y must be vectors, 
t = pv 



t 



net 



r 



3t 



r 






^^pvdV + ^^Q5pv(v • dA) 

J 



net 



or 



(35) 




(36) 



3. Energy 


Equati on 


X = 


pe 


Y = 


Q - IP 




r r 

• 




rr = Q - p 



( 37 ) 
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p 




or 




( 38 ) 



In supplement B, Eq. 38 is transformed into a more 
convenient form by introducing the enthalpy. This form is 
shown below ( read supplement B to see how this was done ). 

r r 

2 



Q - F 



cv 



cv ft 



cs 



( h + 2^^ ^ ^ ^ 



dA) 



(39) 



Where is the net useful (shaft) work. 
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C. SELF-CHECK QUIZ NO. 3 



Select the proper answer 

1. In the generalized notation expression for the FUPLAs 



r 



dt 



XdV 






syst 

The q uanti ty X : 

a. must be a scalar. 

b. is the unknown^of the problem. 

c. represents an extensive property. 

d. none of the above. 



At time "t" the control volume and the system: 

a. are held stationary with respect to an inertial 
frame of reference. 

b. occupies momentarily the same space. 

c. are partially contained one inside the other. 

d. cannot exist at the same time. 

Which of the following expressions is correct based 
on calculus (noting that in the CV, X = X(x,y,z,t)): 






1 im 

At-^-0 At 



1 i m 



1 i m 1 r 
At->0 At 



' 3t 
@t 






CV 



XdV 



/ 




r \ ^ 




syst^^^ - 


syst'^'' - k 


\J 


@t J 


@t+At / 



r 



XdV 



r 

syst 
@ t+At 



XdV 



syst 

@t 



CV 



XdV \= hr 



9t 



XdV 



CV 



lim 1_ 
At->0 At 





f 




d 




. XdV - 
syst 


1 


VJ 

dt 


u 


@t+At 


0t 1 


J 



XdV 



CV 



*(3)C ‘(q)2 ‘(^)l :su0MSUV 
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The following equation relates the expression for 
FUPLAs for a system with an expression valid for 

/V 

the control volume 




a. This equation is directly derived from observa- 
tion of nature. 

b. The first term of the right-hand side is really 

a volume integral and is known as the "convective 
term. " 

c. The second term of the right-hand side of the 
equation is a surface integral sometimes called 
the "convective term." 

d. Can be applied only to unsteady flows. 



• (D )t> :U 3MS u\/ 
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IV . EXTENSION OF JHE CONTROL VOLUME APPROACH TO OTHER 

FUNDAMENTAL PHYSICAL LAWS 



The Law of Angular Momentum and the Second Law of Thermo- 

/s 

dynamics for a system were stated without any further 
considerations in Table III. 

Now we will try to obtain expressions of these laws 

/N /% 

applicable to a control volume. 



A. ANGULAR MOMENTUM 

From the development of the momentum theorem for a system 
in Table II or Eq. 12 






net 



dt 



syst 



V p 



dV 



(40) 



Performing the cross product of both sides of the above 
equations with a vector r whose origin is fixed to the same 
inertial frame of reference relative to which v is measured, 
we can write 



M 



net 



dt 



Jsyst 



f>cp V dV 



(41) 



The quantity in the left-hand side of the above equation 
is the moment of the force with respect to the origin 

of r. We have assigned to the moment the symbol 

The integrand on the right-hand side of Eq. 41 is the 
cross product of the moment arm r and the momentum and defines 
a quantity called moment of momentum or angular momentum. 



Using the generalized notation introduced previously: 



\ 



p(r X v) 



t 



net .. 



(42) 

(43) 
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Figure 10. Coordinate System 



Using the expression for the Fundamental Physical Laws 
developed for a control volume in Eq. 29: 



r 



r 



dt 



J 



syst^f- ' It 






cv 



p ( r X V ) dV+ 



O 



Q 3 p( r X v) ( V dA) 



(44) 



So 







c 


M = ^ 

net 9t 


^yp(r X v)dV +C 


\s ^ ( V‘ dA) 









This is the so-called Moment of Momentum equation. 

B. SECOND LAW OF THERMODYNAMICS 

The Entropy Equation for a system as the mathematical 
expression of the Second Law of Thermodynamics was stated in 
Table III. 
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inequality valid for an irreversible process 



(45) 



T 



> dS 

t 

I equality valid for a reversible process 



Where 6Q is a small quantity of heat transferred during 
an elementary part of the cycle, T is the absolute tempera- 
ture at that point of the boundary, S the entropy, a 



property and s the specific entropy (entropy per unit mass) 









r 






S 




(46) 


Also 












6Q 


= 


area " 


(47) 


• 




j 







where q is the heat flux (amount of heat flowing) per unit 
ti me per unit area 



^ (§A) 

dA ^dt' 



(48) 



Dividing both sides of Eq. 47 by T 



T 



r 

area f 



(49) 



Substituting Eq. 45 into Eq. 49 



dS > J 


area 


q 

T 


dt dA 


(50) 




r 


• 






^ > 
dt - J 


area 


g. 

T 


dA 


(51) 



but from Eq. 46 



r 



dt 



dt ^ 



syst 



sp dV 



(52) 



and finally from Eqs . 51 and 52 



dt 



syst 



s pdV ^ 






area ^ dA 



(53) 
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C. SELF-CHECK QUIZ MO. 4 



1. Can you associate Eq. 53 to the generalized 
notation that we had already developed? 

2. What is the difference, besides the type of 
quantities involved, between Eq. 53 and the 
expressions for the other FUPLAs for a system 
in Table II? 

3. Can you develop an expression for the Second Law 
of Thermodynamics valid for a control volume? 

(See answers on next page) 
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Answers to Self-Check Quiz Mo. 4 



1 . 



X 



PS 

r 



Y > 






9. Hfl 
Area T 



2. The sign ^ depends on what type of process is taking 
place for the heat being transferred in or out of 
the system. The equal sign holds if the process is 
reversible. The inequality sign holds if the 
process is i rre vers i b 1 e . 



3. Applying X and Y to Eq. 29 



7“ 



at 



r 



^yPSdV ^^QsPsC'^ • dA) ^ 



r 



j 



— dA 
area T 
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V. APPLICATION TO PARTICULAR PROBLEMS 



A. CONTINUITY EQUATION 
From Eq . 34 . 



r 



r 



cv 



dV + d(\^ 

at 



dA) 



0 



Applying this equation to the following 

One-dimensional, steady, incompressible 
Fig. 11 . 



case : 

flow as shown in 




Figure 11. Container with unequal inlet and exit areas. 

One-dimensional flow means that the velocity is a function 

of coordinate. There is no change of 

velocity with respect to other coordinates. 

Steady flow means that the velocity is not a function of 



9UI iq 
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'1 






Incompressible flow means that the 
change . 



does not 



Considering the above, we can immediately make the first 
term of Eq. 34 equal to zero. 



Now Eq. 34 becomes 

r 



J 



cs 



p V • dA 



0 



(53) 



p is a constant so we can take it out of the integral in 
Eq. 54. 

Also V is a function of x only so 



->■ 

V 



■r 
V ^ 
X 



And 



B ut 



dA = dA T + dA 3 + dA t. 

X y z 



V ♦ dA = v„dA^ + 0 dA + 0 dA 
XX y z 



So Eq. 54 becomes 

r 






cv 



(54) 



It is only necessary to evaluate the integral Eq. 54 to 
stations 1 and 2 of the control volume since there is no flow 
through the walls 

(p - (p = 0 



Eliminating the subscript x: v = v-, 

Xi I 



and V s Vp 
X 2 ^ 



V, A, 



Vj A^ 



(55) 



In the evaluation of the integral e'quation care must be 
taken with the signs remembering that dA can be positive or 

A 
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negative, depending on the 

(i.e., at the inlet v • dA is a n u a n t i ty ) . 

In the following example we will show the system and con- 
trol volume approach to an unsteady problem. 

B. EXAMPLE OF APPLICATION OF THE SYSTEM AND CONTROL VOLUME 

APPROACH TO THE SOLUTION OF AN UNSTEADY PROBLEM 

A partially evacuated bottle initially containing a mass 
m.| is to be filled with air from the atmosphere until the mass 
reaches a value of Find the work to be done to achieve 

this and if there is no heat interaction, find the change in 
internal energy. 

1 . System Solution 

/N /S /S 

A system or control mass is . 

/s 

We have to choose a system consistent with the above 
definition. The mass during the process must remain unchanged. 
For this purpose we define the system boundaries so as to 
enclose the final mass at all times. (See Figures 12 and 
13). 

We are consistent with the definition because this 
system contains a fixed quantity of mass (m 2 ) although the 
volume that contains it changes. 

The boundaries of the system at time t contain mass m.| 
inside the bottle and mass m 2 - m.| outside the bottle. The mass 

of fixed identity, i.e., the system is 

(write an algebraic expression in terms of m.| and m 2 ). 



*-ui + (*-01 - ^ui) = ^lU 
AqLquopL paxig go sseiu go /Cgigusnb e 

9A iqebau 

eeue 0 L|T buiguesaudeu uogoeA aqg go uon^9“*!-P 
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A 




system 



Figure 12. System solution. State 1. 




i ns ul ated 
bottle 




system 



/s 

Figure 13. System solution. State 2. 
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At time t + At the boundaries of the syste 
changed forming a new volume, but being consistent ,v’th tne 
definition of system, the mass contained in it is still 



Using the First Law of Thermodynamics with the sign 
convention of Fig. 4 we may write 

= Eo - E, '■ (56) 



1^2 



1W2 ^2 ■ 

The first term on the left-hand side of Eq. 56 may be 

eliminated because 

If the only changes in energy are changes in internal 
energy, we can apply the above mentioned expression of the 
First Law of Thermodynamics for a system, reducing it to 



W 

l'^2 



U. 



U 



1 



The change in internal energy is the internal energy 
at t+ At minus the internal energy at t 



AU 



' 2^2 ■ l^'^l^l ^"^2 ■ "^1^ ^atm^ 



(57) 



Where the subscript "atm" refers to the properties 
of the air at atmospheric conditions outside the bottle. The 
atmosphere acts as a reservoir (all properties remain constant) 

The work is the work that has to be done on that 

part of the system outside of the bottle at time "t" to 
introduce it to the bottle. 



W 

1 ^2 



Patm ‘ <"'2 



)] 



(58) 



and Eq. 57 becomes 



Patm^^atm [m, u , - 



(59) 



2 . Control Volume Solution 

If we now solve the same problem using the control 
volume technique, we have to choose the control volume 
conveniently. 
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insulated 

bottlp,^^ 



^atm 




CV — • 

Atmosphere 

Figure 14. Control volume solution 



During the process of filling the bottle the atmos- 
phere is a . Its properties (i.e., ^atm 

) . 



h 



atm‘ 



Recalling Eq. 39, the general expression for the 

Energy Equation for a control volume. 

r 

C) 



CP 



CV 



..-(ep) dV + (h + + z ^)p(v*dA) 

CV'^^ JCS 



In the problem under consideration there is no heat 
transfer and we can neglect both and the potential energy 

changes. Thus Eq. 39 reduces to: 



1 Ir(ep) dV +4^ (h+ -^4jp)p( v*dA) 

JCV -^'CS 'C 



(60) 



We are solving a non-steady flow problem. Integrating 
Eq. 60 with respect to time, t^ being the time at which the 
bottle is filled with mass m^ , we have 



poSueqoun uietuou 

jiOAuasoj ;suoMSUv 
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J 



0 






r 



V 



(ep) dV +(5 (h + 
CV Jcs 



) p { V • dA) 



dt = 3 



( 31 ) 



Working the integral Eq. 6 1 by terms, we have for the 
first term 



rt. 



0 



r 



JZ\ 



. 3T (®P) dV 



dt 






1 



dE 



(62) 



And if all changes in energy inside the Control volume are 
changes in internal energy, Eq. 62 reduces to: 





• 










0 




CV 



dt 



“2 - 



(63) 



Where the subscripts 2 and 1 mean final and initial states 
respectively. Now let us define 



h = h + ^ 

t 2g, 



(64) 



the total or stagnation enthalpy. Replacing Eq. 64 in the 
second term of Eq. 61 we obtain: 



dt 







r 

2 




^^f 




r 


0 


J 


cs * 2 g 


dt = 

J 


0 


J 


h^ p(v-d^) 



(65) 



dm 



But p(vdA) = ^ . For details on this see Fig. 9 and the 
development of Eq. 26. With this, Eq. 65 reduces to: 

^^f 



y 



^ '•^f 




0 





cs 



(h + ^^^)p(v-dA) 



dt 






/, • dm, . dm, 

(ht 3F> - '’dt> 

out in 



dt 

( 66 ) 
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f 






There is no outward flow of mass. The total inflow 
mass during the interval from time 0 to t^ is ( ^ ■ ni-j). 
During this interval, the total or stagnation enthalpy of the 
air entering the control surface is equal to the constant 
atmospheric enthalpy With all these considerations. 

Eg. 66 reduces to: 










r 


T 




(h + 2 )p(vdA) 
CS ^^c 


0 





dt 



- ( m. 



- ) 



atm 



(67) 



Substituting Eqs. 63 and 67 in Eq. 61 we obtain: 



“2 - ^ 



>’atm('"2 - "l> 



And if 



U 2 ^2 



*^atm “ *^atm atm^ "^atm 



U-j = u.| m.| 



Patm 'ra«(™2 - "'1> = - [m,u, - (m^ - m,) 



( 68 ) 



Equation 68 is exactly the same as Eq. 59 obtained using 
the system approach even though a few additional steps were 
needed. 
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VI. SUPPLEMENT A: BODY AND SURFACE FORCES 



In discussing the Momentum Equation we had in Eq. 11 



net 



volume 



I' 



surface 



And at that time the distinction between a body and a surface 
force was not explained in detail. 



The body or volume forces are proportional to either the 

mass or volume of the body and comprise those forces involving 

action at a distance. These forces are a consequence of long 

range force fields, and include such examples as the force of , 

gravitational attraction, magnetic forces and electrical 

(coulomb) forces. For accelerating coordinate systems, "inertia 

forces" such as centrifugal and coriolis forces would also be 

included. But in this case these last two forces must not be 

considered because our development was done for an inertial 

2 

coordinate system. 



Surface forces are those forces which are exerted at the 
control surface by the material outside the control volume on 

/s 

the material inside the control volume. These forces are a 
consequence of short range force fields. Such forces are 
exerted in the form of surface stresses. We can distinguish 
two types of surface forces: (1) those arising from normal 

stresses, or pressures, acting on the control surface, and 
(2) those arising from shear stresses, or viscous stresses, 

/N 

acting on the control surface. 



2 



See footnote on page 22. 






Body or vol ume 
forces 



g ra vi tat i on al 
magnetic 
lelectri cal 
i nerti al 



External 
Forces 
making up 

^net 



Surface forces 



{ pressure -v n ormal 

viscous or ^ 
shear " 



cen tri f uga 1 

Coriolis 



stress 
r stress 
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VII. SUPPLEMENT B; 



INTERMEDIATE STEPS IN THE OERIV"-"!ON 



OF THE ENERGY EQUATION 



Applying the general expression for the FUPLAs. Eq. 31, 
to the Energy Equation 

^ f 

dV +Q) oe(v-dA) = Q - P (38) 



CV 



3( pe) 
9t 






CS 



Also stated vnthout proof was that Eq. 38 can be transformed 
into Eq. 39. To do this, what P is must be considered. In 
Eq. 16, which is the most elementary expression for the First 
Law of Thermodynamics, was stated that P was the rate of work 
done by the System to the surroundings. 

In Supplement A the nature of the forces that may be present 
and dividing them into body and surface forces, was examined. 
Body or volume forces do not appear in P since they do not 
represent work done on the s u rroun di n gs . Rather, they appear 
in the energy term as a potential energy (for most cases of 
interest). We write e as: 



e 



u + 




+ gz + . . . 



Hence, we restrict ourselves to the contribution of the surface 
forces top. These comprise normal and shear stresses and the 
objective here is to single out the contribution of the normal 
stresses which are assumed to be equal to the static pressure 
acting on the fluid. 

Consider an element of fluid of mass pdV entering the 
control volume. 
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■ ' » 






1 - 




*— dX — 



dA 



F.igure B1 . Diagram for calculating flow work 

' work = force x displacement 

= (pdA)dX (Bl) 

rate of work = ^ (pdA dX) = ~ (pdV) (B2) 

Hence dividing by the mass (pdV) we get 
rate of work/lbm = ^ (p/p) = ^ (p*l^) (B3) 

3 

The term p/p is called flow work or flow energy since it 
represents the work to get the fluid into or out of the con- 
trol volume. The important thing to note is that we can 
separate out the contribution of the normal stresses to P by 
simply accounting for the flow work. Since the fluid must 
always perform this work upon entering or leaving the control 
volume, this is a very useful step. 



= 1/p is the specific volume (volume per unit mass, 

3 

i.e., ft /Lb^). Do not confuse this symbol V-with V for 
volume or v for velocity. 
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umTITITITITl'lT 




f* = IP^y + (rate of flow work) (B^) 

In Eq. B4, P^y represents what remains, namely, work clone 
by shear stresses and shaft work. Since friction is usually 
handled indirectly, P^y is commonly labeled as the net useful 
or shaft rate of work. 



Let us see what the above means for a simple compressible 
substance. Suppose a fluid enters the control volume at 
station (1) and exits at station (2) and suppose further that 
all the properties are homogeneous at each station. Now, on 
a per^pound basis, the fluid must do an amount of work P-| 
to enter the control volume, and an amount of work p^ to 
leave the control volume. The work that may be performed is 
still 



l“2 



(-1 

p d'Z^ = 



J 



rz 

1 




so that analogous to Eq. B4, 

^2 

W Q y “ P d - ( P 2 “ p 1 1^1 ) 

J 1 



^2 

^ 1 



t?-dp 




1 

P 



dp 



(this is an important result which is usually written down 
without proof. Show this to yourself starting with a steady 
flow form of the First Law and setting AKE = 0 and APE = 0). 

Now we are ready to arrive at the final form of the First 
Law namely, Eq. 39. Referring back to Fig. 9 we can do this. 
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Figure B2. Two-Dimensional Projection at time t+At 
Equation B2 can be written as 



^ (p dAdx) = ^ (p dA • V At) 



(B7) 



Since for the fluid entering dA opposes v and for the fluid 
exiting, they are in the same direction, we get 

r r 

1 i m 
At-^0 



rate of flow work 



1 _ 

At 



Sin J 



1 _ 
At 



= Q p(vdA) 

CS 

J 
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(B8) 



(B9) 








r 






\ 

i 

! 



Hence , from Eq . B4 



IP 



r 





p ( V • d^) 



And the First Law becomes 

r r 



j 



a(pe) 

cv at 



dv + 



o. 

J 



cs 



[u 





^^]p( v-dA) 
^c 



Now substitute h 



u +^to get Eq. 39. 



(BIO) 



Q - IP 



CV 



(BID 
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VIII. SUPPLEMENT C: EXPRESSIONS FOR THE FUN D," "ENTr L '-'-lYSICAL 

LAWS (FUPLAs) IN A DIFFERENTIAL FORM 



The equations derived in the main text of the handout are 
in integral form. There are many applications where it is more 
useful to begin with a differential form of the equations. 



An 

FUPLAs 



alternate for Eq. 34, 
i s 




3t 



dV +0 

J 



CS 



the general expression for the 



r 

X(v • dA) = y dV 

J cv 



(Cl) 



r 

where Y = y dV 

Jcv 



provided Y is expressible on a per unit volume basis. 

Reference 14 or any other Vector Analysis book explains 
in detail a very useful mathematical tool, the Divergence 
Theorem , 



The Divergence Theorem (also called Theorem of Gauss) 
states that, for any volume in a vector field c , the normal 
component c * n integrated over the area is equal to the 
divergence of this vector field integrated over the volume. 
In mathematical symbols we have: 



r 

O 



\ ‘n dA 



r 



J V 



( V* c ) dV 



Reference 3 states that c niay be generalized to any 

tensor field. Now let us apply the Divergence Theorem to the 
second term of Eq . Cl , 

r C 

(t) X(v • dA) = V-(X v)dV (C3) 






CS 



J cv 
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i 











where it is not difficult to recognize that we have identi-^ied 
C in Eq. C2 with Xv . Note that this quantit,v is a vector if 
X is a scalar or a tensor if X is a vector! 



Replacing Eq. 



J CV 



ft ^ 



C3 in Cl 

r 

+ V • X V d V 

-^CV 



r 



y dv 






(C4) 



Therefore, since the control volume is arbitrary and common 
to all three terms 



31 

9t 



+ 



V-(Xv) 



y 



(C5) 



/\ 

This is the differential form of the Fundamental Physical 

/*s 

Laws using the generalized notation X. We will apply Eq . C5 
to entries in Table II. 



A. CONTINUITY 

X = p and y = 0 



t = 0 

For steady flow 

V*pv = 0 

4 / \ 

and for incompressible flow (p = ). 

Thus 

V • V = 0 



(C6) 



(C7) 



(C8) 



Some books (like Ref. 4, p. 44) goes straight from Eq . 
C6 to C8 without assuming steady flow. ’This can be done if 
the flow is incompressible, since then p does not change 
either with time or position. 
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B. MOME.NTUM 



->- 

X = 


pv 


^ = "^net 


^^net 


per unit volume) 






+ V-(pvv) 


= f . 
net 


(C9) 


Note 


that 


the product vv is 


not defined 


in vector analysis. 



It is a dyadic product (a tensor) and we should be using tensor 
notation. However, V*(pvv) represents a vector [Ref. 15, p. 
200 ]. 



The first term of Eq. C9 can be written as 



3(pv) 

at 




p 



ll 

at 



and the second term can be written 



(CIO) 



V* ( pvv ) 



v(V*pv) + (pv*V)v 



(Cll) 



Replacing Eqs. CIO and Cll in Eq. C9 

V (|^ + V-pv) + (pv-V)v + p|^ = f^g^ (C12) 

The parenthesis of the first term in Eq. C12 is zero from 
continuity Eq. C6 , so Eq. C12 reduces to 







r 3 V , 

* 


(v.v)v] 


^net 




(C13) 


The re 


i s 


a compact 


notation 


often used in 


connection 


to these 


type s 


of 


equations 














D[ ] = 

Dt 


9[ ] + 

at 


(v-V)[ ] 




(C14) 


DLJ. 

Dt 


i s 


called the 


"substantial," "total. 


" or "flow" 


d e r i V a - 


ti ve 


according to 


di f f eren t 


authors. Rewritring Eq. 


C 1 3 using 



the notation introduced in C14 
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which is an equation of motion analogous to m — = ^n°t 
the expression for Newton's second law for a solid, 
such as a cannon ball. A deeper insight is necessary' 
regarding the term on the right-hand side of Eq . C15. This 
term is related to whose meaning was explained in 

supplement A. So we can write 



net 



^volume ^ ^surface 



(C16) 



For the moment we will concentrate our discussion on 



^surface since it is not immediately apparent how to express 
it on a per unit volume basis. The surface forces depend on 
the rate at which the fluid is strained by the velocity field 
The system of forces determine a state of stress and the re- 
lationship between stress and strain can only be given 
empirically. Recalling again Supplement A, a distinction 
can be made between forces arising from normal stresses and 
forces arising from shear stresses. 



^surface ^normal ^ ^shear 

The force due to the normal stress may be written 



(C17) 



^normal ~ - V p 

and the force due to the shear stress, under the Stokes 
hypothesi s , 



(C18) 



I,, ear = w v (C19) 

Where p is the viscosity of the fluid. It is a property of 
the fluid (and may be strongly dependent on temperature). 
Replacing Eqs. C19 and C18 in C17 



- 5 - 

f. 



- V p 



p V 



(C20) 



surface 

We have only shown between Eqs. C16 and C20 how the different 

-y- 

kinds of stresses make up ^ complete derivation 

of Eq. C20 can be found in Ref. 4. In the derivation of Eq. 
C20 it was assumed that the fluid is isotropic, namely, the 
components of the stresses are equal in all direction and 
that it is Newtonian, i.e., a linear relation exists between 
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i 



the rate of strain and the stress, the relation bet'/c-en tb" 
two being emp i r i ca 1 . 



Replacing Eg. C20 in C16 




f 



volume 



V p + V 



(C21) 



and replacing Eg. C21 in C15 




f 



VO 1 ume 



2 

V p + y V V 



(C22) 



Eguation C22 is a very useful form of the eguations of motion 
called the Navier-Stokes Eguations (a set of three eguations, 
one for each coordinate). They form the starting point for a 
course in boundary layer theory. Written as in Eg. C22 it 
has only the restrictions of isotropic, incompressible, New- 
tonian flows. The incompressible restriction arises from the 
fact that we had treated y as a constant taking it out of the 
Laplacian in the second term of the right-hand side of Eg. 

C22. For compressible flow, where there exists gradients of 
temperature that affect strongly the value ofp, the above 
mentioned term should be written V*yVv. In the derivation of 
Eg. C22, the compact form of Navier-Stokes eguations, most of 
the effort was given to explaining the significance of each 
term and its origin. For a complete and rigorous derivation 
of these eguations we suggest the use of Ref. 4 in the con- 
text of an appropriate course. 

In Aeronautical Engineering, the Navier-Stokes eguations 
are useful in flows (not fluids) where the viscous effects are 
important, i.e., flows where large velocity gradient exist. 
Such situations are found in boundary layers where there is a 
change of velocity between a value of zero at the wall and a 
value different from zero at the free stream. The gradient 
creates important viscous forces taken i'nto account in the 
last term of Eg. C22. In the free stream outside the boundary 
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layer, although j(the viscosity) exists as in the boun^^'^y 
layer region (same fluid), the viscous forces do not irist 
because there is no velocity gradient, and the last term in 
Eq. C22 is zero. 

In which other case would you do the same ? 

In these two cases Eq. C22 reduces to 

4 ' ^volume - ’P 

Write Euler's Equation of Motion for one dimensional flow 
using the definition of total derivative 



(C24) 



For steady flow and neglecting the volume forces, Eq. 
C24 reduces to: 




1 Ie. 

p dx 



(C25) 



Why does Eq. C25 only contain total derivatives? 



Now integrating Eq. C25 we get 
.2 



V 

2 



r 

^ = 
p 



Or if the flow is 
be integrated to give 
..2 



^ + 



e. 

p 



constant 



constant 



(C26) 



the pressure term can 



(C27) 



(queqsuoD sl d) a[qLss 0 udujoouL 

•MO[g A'peaqs [Buoisuaiuip auQ ‘(aLduiexa 
joj. X /Ces) aqeuLpuooD auo go suoLqoung aue d pue a asneoaq 



Jip-X [OA 



xe 



:! •+ -Je = ( 



,xe 

Ae' 



+ — )d 
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This is called Bernoulli's Equation. Sometimes it has an 
g 

additional term coming from the body forces '.hen v;e 

neglected in Eq. C24 (as they pertain to a gravity field). 

This equation looks like an Energy Equation but it is 
not complete as we shall see in part C of this Supplement. 

Exe rci se : 

What does it mean when the flow is Newtonian? 

On what property does the viscosity primarily depend? 

It can be shown that the assumption of one-dimensional 
flow can be replaced by the restriction that the constant in 
Eq. C27 takes different values for different streamlines. 

In other words, the left-hand side of Eq. C27 is equal to a 
different constant, depending on which streamline of the flow 
is considered. 

Also, it can be shown that Eq. C27 is valid everywhere 
(no matter if we change streamline) if the flow is 

irrotational (Vxv=0). 

0 

Equation C27 is known as Equation. The 

assumptions and constraints in its derivation are: 
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C. ENERGY 



In our development of the Energy Equation for a control 
volume, Eq. 38 we stated that this form was not very con- 
venient and in Supplement B we transformed it into Eq . 39. 

The differential form of the Energy Equation, Eqs. 38 
or 39 can be derived in a manner similar to the Momentum 
Equation, thus only the results will be given here. For a 
description of the various terms and their application see 
Ref. 5. 



The following form is usually quoted [Ref. 5] for the 
Energy Equation: 

*^Ft ■ PV'V = -V*q + p<I> (C29) 

or in terms of the enthalpy 

p|y + p(v*V)h + [||- + (v‘V)p] = -V*q + (C30) 



where e 

Q 





J cs 



(neglecting changes in kinetic 
and potential energy) 



(C31) 



and $ = dissipation function acting through the fluid 
vi scos i ty 



If neglecting radiation heat transfer and assuming that 
Fourier's Law of Heat Conduction applies, then 

->■ 

q = - k V T (C32) 

and V-q = - V(k V T) (C33) 

4 o 

V*q = - k T for k = constant (C34) 

The solution to the equation above represents the effect of 
fluid motion on heat transfer and is studied under the 
title of "Convective Heat Transfer." 
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Substituting Eq 
p|y + p(v-V)h 



C34 into Eq. 31 

+ (v*V) p] = kV^T 



( C35) 
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D. SELF-CHECK QUIZ NO. 5 



1. As a consequence of steady flow 

a. the velocity is not a function of the position of 
the f 1 ui d parti cle . 

b. density is not a function of the velocity. 

c. velocity is not a function of time. 

d. density is constant. 

2. What term(s) would you eliminate from Eq. C22 if you 
are told that the fluid/ flow is 

a . I n V i s c i d . 

b . Incomp ress i b 1 e . 

3. If for a fluid element it is found that there exists 

a linear relation between the rate of strain and stress, 
the fluid is said to be 

a . E u 1 e r i a n . 

b . Newtonian. 

c. Lagrangian. 

. d. Isotropic. 

4. From the differential form of the Momentum Equation we 
derived several equations associated with great 
scientists and mathematicians of the past. The expres- 
sions known as Euler, Bernoulli and Navier-Stokes 
Equations have in common that they are: 

a. Energy Equations. 

b. Continuity Equations. 

c . Moment Equati ons . 

d. Equations of Motion. 
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5. Match with the name of the equations (Column 1) (ior-'v^d 
in this supplement all the assumptions or constraints 
made in its derivation (Column 2). 





Column 1 




Col umn 2 


a . 


Navier-Stokes Equations 


a . 


I n V i s c i d fluid. 


b. 


Euler Equation 


b . 


One-dimensional flow. 


c . 


Bernoulli Equation 


c . 


I sot ropi c fluid. 






d. 


Incompressible flow. 






e . 


Newtonian fluid. 






f . 


Steady flow. 






g- 


No velocity gradients. 



‘ ^ ‘ p ‘ q * 6 uo e 

• 6 JO B 
• 0 ‘ p ‘3 
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